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Lifetime measurement of the 6792 keV state in 15 important for the astrophysical 

S factor extrapolation in 14 N(p,7) 15 

D. Schiirmann/'BR. Kunz, 1 I. Lingner, 1 ^ C. Rolfs, 1 F. Schumann, 1 '! F. Strieder, 1 '! and H.-P. Trautvetter 1 

1 Institut fur Experimentalphysik III, Ruhr-Universitat Bochum, Germany 
(Dated: April 11, 2008) 

We report on a new lifetime measurement of the E x — 6792 keV state in 5 via the Doppler-shift 
attenuation method at the E = 259 keV resonance in the reaction 1 N(p,7) O. This subthreshold 
state is of particular importance for the determination of the ground state astrophysical S factor 
of 14 N(p,7) 15 at stellar energies. The measurement technique has been significantly improved 
compared with previous work. The conclusion of a finite lifetime drawn there cannot be confirmed 
with the present data. In addition, the lifetime of the two states at E x = 5181 and 6172 keV have 
been measured with the same technique in order to verify the experimental method. We observe an 
attenuation factor F(r) > 0.98 for the E x = 6172 and 6792 keV states, respectively, corresponding 
to r < 0.77 fs. The attenuation factor for the E x = 5181 keV state results in F(r) = 0.78 ± 0.02 
corresponding to t — 8.4 ± 1.0 fs in excellent agreement with literature. 

PACS numbers: 26.20.Cd, 24.30.-v, 26.65.+t , 27.20.+n 

Keywords: lifetime measurement, hydrogen burning, Dopplcr-shift method, gamma-spectroscopy 
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I. INTRODUCTION 



The capture reaction 14 N(p, 7 ) 15 (Q = 7297 keV, fig. 
[1]) is the slowest process in the hydrogen burning CNO 
cycle [l| and thus of high astrophysical interest. This 
reaction plays a key role for the energy production of 
more massive main sequence stars and the detailed un- 
derstanding of the neutrino spectrum of our sun [2j, |3| 
as well as the age determination of globular cluster stars 
[4(. The reaction was recently studied in three exper- 
iments at energies ranging from E — 70 to 480 keV 
[a, |g, 0, @] (all energies are given in the center-of-mass 
frame if not indicated different) and previously over a 
wide range of energies, i.e. E = 240 to 3300 keV (0 
and references therein). The major motivation for rein- 
vestigating this reaction was a reanalysis [lfJ| of the data 
of [9j using R-matrix resulting in a significant change in 
the extrapolation of the ground state astrophysical S fac- 
tor from Sg S (0) = 1.65 keV b to a nearly negligible 
contribution. A reduction of the ground state contribu- 
tion could indeed be verified by the recent experiments 
[a Q. However, the R-matrix analysis revealed that be- 
low the Er = 259 keV J 71 = l/2 + resonance the data 
followed primarily the low energy wing of the resonance 
(fig. [2]) and could not probe significantly the behavior 
of the interference structure of the four J 71 ' = 3/2 + reso- 
nances at En = —507, 987, 2187 keV and the background 
pole - here fixed at 6 MeV - around the energy of this 
J 71 " = l/2 + resonance which in turn adds incoherently to 
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FIG. 1: Level structure of 15 



the ground state S factor (for details see [a] and PJJ). A 
precise knowledge of this interference structure is needed 
for reliable extrapolation due to a minimum of the S fac- 
tor curve near E = 160 keV. 

In a second experiment at the LUNA accelerator at 
the Laboratori Nazionali del Gran Sasso, Italy, with a 47r 
BGO detector and a windowless gas target system [7| the 
total S factor of 14 N(p,7) 15 has been measured down to 
E = 70 keV with 5(70) = 1.74 ± 0.14 keV b (statistical 
error only). However, the ground state contribution at 
that energy is expected to be S gs (70) = 0.07 keV b and 
falls rapidly above 70 keV up to the resonance [f|. It is 
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FIG. 2: Astrophysical S factor of 14 N(p, 7 ) 15 0. Shown are the 
measured total S factor (full circles, [3]), the ground state 
S factor of Formicola et al. (open circles, |5|]) and Runkle et 
al. (stars, |8J]), and the R-matrix calculation of 5] for the total 
S factor (solid line) and the ground state transition (dashed 
line) as well as the ground state R-matrix calculation of [12i ] 
(dotted line). 



obvious that the experiment [7J cannot probe sensitively 
the contribution of the ground state transition within 
a factor of two of the results of J5j. Fig. [2] shows a 
comparison of the available data together with R-matrix 
fits [a, [||]- Note: the data of [7| are not taken into account 
in the R-matrix analysis. 

The energy dependence of S ffs at low energies is pri- 
marily controlled by the width of the subthreshold state 
at E x = 6792 keV. The width could be deduced by a mea- 
surement of the lifetime of this state. In this case such an 
experiment is feasible by means of 7 spectroscopy if the 
lifetime is long enough to allow for a observable Doppler- 
shift attenuation in the slowing down process of the re- 
coiling 15 nucleus. This method is called Dopplcr-Shift 
Attenuation Method (DSAM) and was reported in [12| 
with an attenuation factor F(t) = 0.93±0.03 correspond- 
ing to a lifetime of r = 1.6 ± 0.8 fs. In this experiment 
the authors used three Ge detectors placed at 0°, 90°, 
and 144° with respect to the beam axis. The advantage 
of this approach is that all information in principal could 
be obtained in a single run, but results in the difficulty 
that each spectrum had to be calibrated absolutely and 
as a function of time to a high precision (better than 
±0.5 keV) in order to extract an F(t) value near unity. 

In contrast to the previous experiment 12] the present 
approach uses one detector placed on a turntable with 
an angle span from -40° to +116°. Gamma spectra were 
obtained at 11 angles of which 2 angle settings were sym- 
metric around the beam axis leaving 9 different angles. 
In this way a left/right asymmetry in the detector place- 
ment could be detected as well as a possible asymmetry 
in the beam alignment. Therefore systematic uncertain- 



ties could be greatly reduced. The gain stability of the 
electronics in the course of the experiment was checked 
by the E 1 = 6129 keV 7-ray emission from the reac- 
tion 19 F(p,a7) 16 0. Due to the fact that all measurements 
were done with the same detector a precise absolute de- 
termination of the 7-ray energy was not necessary. The 
determination of the slope of the energy calibration was 
obtained via the 14 N(p,7) 15 reaction itself using the 7- 
ray energy information of the various transitions in 15 
given in [6J. 



II. EXPERIMENTAL SET-UP AND 
EXPERIMENTAL PROCEDURE 

The 500 kV accelerator at the Dynamitron tandem lab- 
oratory at the Ruhr-Universitat Bochum provided the 
proton beam with typically 100 to 150 [iA on target at 
E p = 318 keV. The energy calibration and the energy 
spread of the proton beam are known to better than 
1 keV, and equal to 70 eV, respectively [la ). Fig. [3] 
shows a schematic diagram of the experimental set-up. 
The beam passed through two Ta collimators with a di- 
ameter of 10 and 4 mm at a distance to the target of 80 
and 6.5 cm, respectively, and was stopped in the N tar- 
get. A LN2 cooled shroud placed at 20 cm in front of the 
target avoided C-build up and ensured a typical pressure 
in the target chamber of about 5 x 10~ 7 mbar. After 
focusing the beam, the 4 mm collimator was removed to 
minimize beam induced background originating from the 
collimator. The target chamber was formed by a cylin- 
drical T-shape steel tube (inner diameter = 9.7 cm, wall- 
thickness = 2 mm) in order to ensure the same influence 
of the 7-ray absorption in the chamber walls for all an- 
gles. The target was mounted on a directly water cooled, 
rotatable target holder. On the target holder three tar- 
gets could be placed at the same time and brought into 
the beam without breaking the vacuum. The alignment 
of the target chamber and the target itself was checked 
to be better than ±1 mm. 





FIG. 3: Schematic diagram of the experimental set-up. In 
the upper left corner the directions of the geometrical offset 
parameter a and b of the turntable axis are indicated (for 
details see text). 
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FIG. 4: Thick-target 7-ray yield curve of the E R — 259 keV 
resonance (proton energy E p — 277 keV) of N(p,7) 1 O ob- 
tained with a 14 N implanted Ta sheet (220 keV 14 N ion beam). 
The proton energy E v — 318 keV for all the runs is indicated 
by the arrow. 



A turntable was placed at the target chamber with 
its axis at the target position. The 0° position could 
be aligned to the beam axis within ±0.5°. The total 
uncertainty of the turntable axis position in respect to the 
beam spot position could be estimated by ±3 mm, where 
the uncertainty coming from the detector positioning on 
the turntable arm as well as the beam spot position on 
the target are included. The angles reached by this set- 
up ranged from -40° to 116° with respect to the beam 
axis (fig. [3]) and reproducibility within less than 0.2° 
was varified. A HPGe gamma detector was placed on the 
arm of the turntable (108% relative efficiency, resolution 
= 2.1 keV at E 7 — 1.3 MeV, front face distance to the 
target = 13.1 cm). 

The N target has been produced by implanting 14 N 
into a Ta backing with a 0.5 mm thickness (backing di- 
ameter = 50 mm). A molecular N beam at an energy 
of Ei4 N2 = 220 keV from the 500 kV accelerator was 
rastcred with a pair of magnetic x-y scanners over an 
aperature (diameter = 24 mm) in close geometry to the 
target leading to an implanted area of 25 mm diameter. 
The 14 N distribution in the Ta was investigated using the 
En = 259 keV resonance of 14 N(p,7) 15 0. A target pro- 
file obtained after an 14 N irradiation dose of 3.9 x 10 18 
atoms/cm 2 is shown in fig. [4] As seen in the figure the 
14 N atoms are distributed from the surface of the Ta to 
a depth of 46 keV at half maximum. A proton energy 
of E p — 318 keV, where the thick-target yield plateau 
is still almost constant, was chosen for all the measure- 
ments. The target stoichiometry after the implantation 
reaches saturation should be equal to that of the com- 
pound Ta 2 N 3 Hi El . 

In addition, a CaF target has been evaporated on a Ta 
backing (thickness of the CaF layer = 410 nm) and was 
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FIG. 5: (a) Gamma-ray spectrum near the E^ — 6129 keV 
transition curve of 19 F(p,cry) 16 obtained with the CaF tar- 
get at E p = 318 keV. The solid and dotted lines represent 
measurements at 0° and 90°, respectively, and demonstrate 
the effect of the Doppler-shift at different angles. The inset 
(b) shows the stability of this 7-ray line versus time during 
the experiment in respect to the average value for both detec- 
tion angles. The data are corrected for the time dependent 
energy shift. 



used for calibration purpose. This CaF target was placed 
on the target holder together with the 14 N target and 
moved into the beam with the proton beam untouched 
in both focus and energy. This procdure was repeated be- 
fore and after each run on the 14 N target. Thus the gain 
stability of the detector/electronic system was controlled 
by observing the E 1 = 6129 keV line of the 19 F(p,a7) 16 
reaction at E p = 318 keV (E = 302 keV in the 19 F + p 
system) at 0° and 90°, respectively. The peak position 
of the Ej — 6129 keV line in each calibration run was 
extracted from the spectrum by the determination of the 
centroid of the peak area. Sample spectra are shown in 
fig- EI The effect due to different contributions from un- 
stopped 16 nuclei is clearly visible. The absence of a 
high energy tail suggests a strongly foreward peaked an- 
gular distribution of the a-particles and as a consequence 
one observes at 0° besides the 7-ray emission from 16 
nuclei at rest a low energy tail due to the backward mo- 
tion of the 16 nuclei (solid line in fig. [5J. However, at 
90° also a high energy tail emerges. The treatment of the 
background below the peak from the unshifted contribu- 
tions depends on the kinematics and angular distribution 
of the a-particles which is not precisely known. Differ- 
ent treatments of this background lead to different cen- 
troids which can amount up to 1 keV difference between 
a 0° and a 90° measurement. We have adopted a linear 
background for the 90° run and an error function (po- 
sition and width taken from the peak) for the 0° case. 
The centroid for both runs agreed within 0.15 keV us- 
ing these assumptions. The result of the calibration runs 



is that the offset in the energy calibration was constant 
within ±0.08 keV during the course of the experiment 
(about 10 hours) while the relative error in the slope was 
5 x 10~ 5 . Moreover, the slope shows a weak linear time 
dependence of about 0.025 keV/hour. All spectra were 
corrected for this effect. In summary, the accuracy of the 
energy calibration is of the order of ±0.15 keV (see inset 
in fig. ED. 

Finally, the angles measured with the N target were 
chosen to be 90°, 0°, 105°, 75°, 20°, -20°, 116°, -40°, 
60°, 40°, 82°, 116°, 0°, and 90° in this order in time. The 
target holder could be placed in two positions: at +110° 
(1) and at -110° (2) with respect to the beam axis. The 
first 11 angles were obtained with the target in position 
1 while the last three angles were measured in position 
2. The different target positions were used to exclude 
possible influences of 7-ray absorption through the target 
holder on the determination of the energy-centroid of the 
peak. Moreover, the fine grid of angles allows to check 
experimentally a displacement of the turntable in respect 
to the beam spot position. 

In each run or angle measurement, respectively, the 
centroid of five 7-ray lines was determined: the sec- 
ondary transitions of the E x = 6792 keV state, primary 
(E y = 1384 keV) and secondary of the E x = 6172 keV 
state, and primary (E^ = 2385 keV) and secondary of 
the E x = 5181 keV state. The primary transition into 
the E x — 6792 keV state was not taken into account 
due to too small shifts of this line. The observed 7-ray 
energy, Ey(0), at an angle 9 in respect to the recoil di- 
rection is related to the emitted 7-ray energy , E° (in the 
center-of-mass frame), by the relation: 



vo 



E 7 (9) =E°-(1 + ^Lp.F{ T )cos{e)) 



(1) 



where ^o is the initial recoil velocity and P is a correction 
accounting for the finite size of the Ge detector, which 
has been obtained by a GEANT4 simulation 16] to P = 
0.987. 

The quality of the spectra of the present experiment 
arc displayed in fig. [6] for the most important transi- 
tion: the E x = 6792 keV state secondary transition. The 
statistics of the 90° run (filled circles in fig. [6]) was im- 
proved by a factor two in the present experiment com- 
pared to previous experiments, i.e. [l2|. The background 
of the E~f — 6792 keV 7-ray line was - similar to the 
calibration runs - estimated by an error function and is 
nearly negligible for this line. We have not analysed the 
single escape lines - which in principal could be also used 
- because they have less statistics and higher background, 
which would in turn introduce larger systematic uncer- 
tainty due to the background treatment. The width of 
the region for the centroid determination was kept con- 
stant for all angles but shifted with the Doppler-shift of 
the line. 

The analysis of the E^ = 5181 keV line was seriously 
hampered by the double escape line (2? 7 = 5150 keV) 
of the 6172 keV peak. Both lines shift differently and 
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FIG. 6: Gamma-ray spectra near the E 1 = 6792 keV tran- 
sition curve of 14 N(p,7) ls O obtained at E p = 318 keV. The 
stars and the filled circles represent measurements at 0° and 
90°, respectively. The lines are to guide the eye. 
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FIG. 7: Gamma-ray energy of the Doppler shifted E-y = 
6792 keV line versus the detection angle 9. The data points 
(filled circles) and the best fit (F(t) = 0.99, solid line) of 
the present experiment are compared with the results (open 
squares) and the fit (F(t) = 0.93, dotted line) of Bertone et 
al. Q. 



therefore at 0° (116°) the distance of both lines was only 
29 keV (34 keV) complicating the background determi- 
nation. These circumstances introduce additional sys- 
tematic errors in the determination of the centroid for 
this transition. The intensity of the primary transition 
7-ray lines is much higher and hence the precision with 
which their centroid could be analysed relative to the ex- 
pected full shift is nearly the same as for the high energy 
lines. The positions of these lines, which are expected 



to show full Doppler-shift, serve as a constraint for the 
determination of the geometrical offset parameters of the 
turntable relative to the target position. 

The experimental shifts A are obtained from the slope 
of a linear fit to the observed 7-ray energies versus cos(6), 
with 8 being the observation angle. The attenuation 
factor F(t) = A/(Af u u-P) can be calculated analyti- 
cally from the measured shift values A. A correction of 
the observation angle 9 due to the real position of the 
turntable was applied, where the geometrical shift pa- 
rameters a and b were kept as free parameters in the fit. 
In a combined fit of all lines using eq. 1 the F(t) val- 
ues as well as the parameters a and b were determined, 
where the F(r) values for the primary transitions were 
only allowed to vary within ±0.01 from unity. The best 
fit results in an attenuation factor of F(r) = 0.99 ± 0.01 
for the E 1 = 6792 keV 7-ray line and is displayed in 
hg- El (solid line) together with the experimental data 
points (filled circles) for this transition. For comparison 
the measured values (open squares) of [12| and their best 
fit with F(t) — 0.93 (dotted line) are shown: the mea- 
surements are in good agreement with the present data 
but the attenuation value of [12j is excluded by more than 
5cr by the dense experimental data grid of the present ex- 
periment. The conclusion from the present experiment 
is that the result is in good agreement with full Doppler- 
shift. Therefore, only a lower limit of the attenuation 
factor F(t) for the E x — 6792 keV transition can be 
given and, in turn, only an upper limt for the lifetime 
of the associated state (see below). The results for all 
transitions are summarized in table U together with the 
calculated uncertainties and the reduced x 2 - The resid- 
uals of all 5 fits are shown in fig. [8] The geometrical 
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FIG. 8: Residuals of all analysed 7-ray lines for the best com- 
bined fit of the Doppler-shift attenuation according to eq. 1. 



offset of the turntable was determined to a = —2.6 mm 
and b — 0.4 mm within expectations. 



III. LIFETIME DETERMINATION 

For the lifetime extraction we followed the assumption 
of PJi, [l8| that the implanted N ions occupy interstitial, 
rather than substitutional, sites within the Ta lattice and, 
therefore, the density should be close to that of Ta. In 
this case F(t) can be calculated by 



F(t) = 



f °° e-*/^(t)cos(^)di 



VqT 



(2) 



We have run the SRIM program [19j to obtain the 
collision details for a large number of recoiling 15 nu- 
clei from which v(t)cos(ip), the time dependent averaged 
projection of the recoil velocity distribution on the beam 
axis, was extracted and then fitted by a sixth order poly- 
nomial. Eq. 2 can then be solved analytically. The 
resulting F(r)-curve is not as steep as found in [I2| (see 
also open symbols compared to solid line in fig. [5]) but 
results in a lifetime determination of the E x = 5181 keV 
state in 15 (table Hand solid point in fig. [5]) which is 
in excellent agreement with Jljj (vertical lines in fig. \§§. 
The analysis leads to an upper limit of the lifetime of 
r < 0.77 fs for the E x = 6172 and 6792 keV states, indi- 
cated by the shaded area in the upper left corner of fig. 
|51 For comparison, fig. [5] also shows the values obtained 
by [l2| as open symbols. If the saturation of N in the 
Ta was not reached completely during the implantation 
the resulting F , (r)-curve is not as steep as in the pre- 
sented case. Consequently, in the worst case of a N:Ta 
ratio of 1:1 the \-a upper limit increases by only 2 % for 
the E x = 6172 and 6792 keV states and to r = 9 fs for 
the E x = 5181 MeV state within the given uncertainty 
range. 

In summary, a finite lifetime for the latter states in 15 
cannot be extracted from the present data. This leads to 
the final conclusion that there are strong doubts that 
the lifetime of these states can be measured by means 
of the Doppler-shift attenuation method. Moreover, ref. 
[18| opens the possibility of an independent check of our 
analysis procedure on experimental F(t) values. In this 
work the lifetime of states in 14 N were obtained using 
the Doppler-shift attenuation method with the reaction 
13 C(p,7) 14 N at E p = 1150 keV. The authors also evalu- 
ated the lifetime with Monte Carlo calculations and in- 
vestigated in the course of the measurements experimen- 
tally the stopping power of 15 N ions in Ta. The reason to 
choose 15 N instead of 14 N was only due to the absence of 
any narrow 14 N resonances in the range of their acceler- 
ator. Our analysis using the collision details (see above) 
obtained from SRIM [19| is in very good agreement with 
the results of [18| , see table |TTJ We also have reproduced 
the experimental depth profiles of 15 N ions implanted in 
Ta of [18] using SRIM pj|, indicating that the stopping 



TABLE I: Experimental F(t) values and corresponding lifetimes compared with literature. 



E 1 [keV] 



F(t) 



X 



r[fe] 



F{rY 



t [is] a 



[is] 



1384 c 
2375 c 
5181 d 
6172 d 
6792 d 



1.00 ±0.01 
0.99 ±0.01 
0.78 ±0.02 
0.99 ±0.01 
0.99 ±0.01 



0.6 
1.0 
0.4 
1.1 

0.7 



3.4 ±1.0 

< 0.77 

< 0.77 



0.68 ± 0.03 
0.91 ±0.05 
0.93 ± 0.03 



9.67 ±1.30 
2.1 ±1.3 
1.6 ±0.8 



8.2 ±1.0 
< 2.5 
<28 



°ref. hj 
6 ref. [Hf 

c primary transition 
^secondary transition 
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FIG. 9: F(t) function according to eq. 2 for ls O in Ta. 
The allowed region and the results of the present experiment 
are indicated by the shaded area in the upper left corner 
for E x = 6172 and 6792 keV and the filled circle for the 
E x = 5181 keV transition, respectively. The results of [13 ] 
are shown as open symbols (circle = 5181 keV, triangle = 
6172 keV, and diamond = 6792 keV transition). The lit- 
erature value for the E x — 5181 keV state lifetime [17J ] is 
indicated by the vertical lines. 



TABLE II: Lifetime determination with the present analysis 
method from experimental F(t) values [lQ] obtained by the 
reaction 13 C(p,7) 14 N compared with the results from [18| . 
The very good agreement of the results demonstrate the va- 
lidity of the present analysis for the lifetime r. 



level in 14 N 
E x [keV] 



ref. [18] Bister et al. present 

F(t) t [fs] t [is] 



3948 
5690 



0.894 ± 0.005 
0.82 ± 0.08 



5.4 ±0.4 

16 ±8 



5.1 ±0.4 
14 ±8 



power calculation from [19( is sufficiently precise for the 
present experiment. 

IV. SUMMARY AND CONCLUSION 

The present experiment has improved the statistical as 
well as systematic uncertainties in the determination of 



the Doppler-shift attenuation of three 7-ray transitions 
in 15 0: 5181, 6172 and 6792 keV. Our upper limits on 
the lifetimes of the E x = 6172 and 6792 keV is only 
marginaly in agreement with the conclusion of [12j , while 
the experimental data of [12| are not in contradiction 
with the present data. 

The upper limit for the lifetime of the E x — 6792 keV 
subthreshold state in 14 N(p,7) 15 converts to a lower 
limit of the 7- width of this state of T 7 > 0.85 eV. This re- 
sult is in good agreement with the value found in [5j with 
T 7 = 0.8 ± 0.4 eV from R-matrix analysis of the S factor- 
curve of the ground state transition in 14 N(p,7) 15 0. The 
S factor extrapolation from the R-matrix fit of [8| indi- 



cates a value of r. 



1.6 eV not in contradiction with 



the present result, but in disagreement with [5j (Note 
the R-matrix fit of [8J] was obtained by omitting the high 
energy data of Schroder et al. [9|]). In a Coulomb exci- 
tation experiment [2fJ| a width of T 7 = 0.95^0 95 e V was 
reported in agreement with the present result. Hence, the 
statements in [a, |6|, |8(] concerning the ground state tran- 
sition in 14 N(p,7) 15 are essentially unchanged, but fur- 
ther conclusions need to await a more detailed R-matrix 
study. Moreover, further investigations for the ground 
state S factor of the reaction 14 N(p,7) 15 are still neces- 
sary. Such an experiment [111 ] was performed recently, i.e. 
a precision measurement of the ground state transition 
above the Er = 259 keV resonance with a Clover detec- 
tor, reducing the necessary summing correction by one 
order of magnitude compared to previous direct studies 
[3, Q • The results [21| indicate again a very good agree- 
ment with the present result and the T 7 of Formicola et 
al. [|. 
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